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The density matrix renormalization group method for ladders works much more efficiently with 
open boundary conditions. One consequence of these boundary conditions is groundstate charge 
density oscillations that often appear to be nearly constant in magnitude or to decay only slightly 
away from the boundaries. We analyse these using bosonization techniques, relating their detailed 
form to the correlation exponent and distinguishing boundary induced generalized Friedel oscillations 
from true charge density waves. We also discuss a different approach to extracting the correlation 
exponent from the finite size spectrum which uses exclusively open boundary conditions and can 
therefore take advantage of data for much larger system sizes. A general discussion of the Friedel 
oscillation wave-vectors is given, and a convenient Fourier transform technique is used to determine 
it. DMRG results are analysed on Hubbard and t — J chains and 2 leg t — J ladders. We present 
evidence for the existence of a long-ranged charge density wave state in the t ~ J ladder at a filling 
of n = 0.75 and near J/t « 0.25. 



I. INTRODUCTION 

Bosonization analyseJli'i'i tcgether with finite sizg 
numerical work using the DMRGdu and other methodsQ 
have given a clear understanding of the behavior of the 
2-leg Hubbard and t — J ladder models. The Hubbard 
Hamiltonian is written: 

^ = E + h.c.) (1.1) 

~ ^E (42aC»lQ + h.c.'j + U ^n,,A,T«»,A,i 

Here c^aq destroys an electron on rung i and leg A = 1, 2 
with spin a =t, i- rii^x^^ is the electron number operator. 
The t — J Hamiltonian is: 

H= - (cI+iAaC»Aa + h.c}j 

- *E (c^2aC^la + /i.e.) 
i.a 

+ J ^ (^SiX ■ Si+i\ 

I, A 

+ JE(5^-^--™) (1-2) 

i 

where Six = cjA ^/^ '^i>' ^'^'^ Hilbert space now ex- 
cludes all states with doubly-occupied sites. We will 
generally set t = 1 in what follows. Both models are 
expected to be in a "CISO" phase, over a wide range 
of parameters, in which the low energy degrees of free- 
dom consist of a single free massless charge boson, whose 



excitations carry even multiples of the electron charge. 
This is often characterized as a "d-wave superconductor" 
based on the nature of the power law decay of pair cor- 
relations which have a positive sign for singlet rung-rung 
or leg-leg correlations and a negative sign for rung-leg 
correlations. Both bosonization and DMRG results exist 
on multi-leg ladders with analytical and numerical un- 
certainty which increase with the number of legs. 

Most of the DMRG data is obtained with open bound- 
ary conditions along the legs of the ladder. Unlike with 
periodic boundary conditions, this generally leads to 
charge density oscillations in the groundstate. If these 
oscillations persist at the center of the chain for arbitrar- 
ily long chain length then they correspond to a charge 
density wave (CDW). More commonly, they decay away 
from the boundaries with a power law, in the limit of an 
infinite chain. In this case, we may think of them as gen- 
eralized Friedel oscillations where the chain ends them- 
selves act as impurities which induce gradually decaying 
density oscillations. Unlike in a Fermi liquid, these gener- 
alized Friedel oscillations in a Luttinger liquid decay with 
an exponent which depends on the interaction strength, 
and which, js simply related to the density correlation 
exponentE'EI Furthermore, the wave-vector of the oscil- 
lations itself can be changed by the interactions. These 
density oscillations, for the largest systems studied, up to 
eight legs, have been identified with "stripes" S3 An un- 
derstanding of the occurance of stripes in these models, 
whether or not they require long range Coulomb inter- 
actions to exist and their connection with superconduci- 
tivity are important open problems. Since any ladder 
system is ultimately one-dimensional if the number of 
legs is held fixed and the length taken to oo, the Friedel 
oscillations /CDW analysis may provide the appropriate 
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description of stripe behavior, at least in this hmit and 
in cases where the stripe wave-vector is parallel to the 
chains. 

In Sec. II we briefly review the bosonization picture 
of the 2-leg t — J ladder. We also present a technique 
for extracting the correlation exponent purely from data 
with open boundary conditions which has not, to our 
knowledge, been used previously on ladders, although it 
has been used on chains. This method is used to obtain 
values for the charge velocity and exponent. In Sec. Ill 
we review the bosonization treatment of commensurate 
charge density waves and the related behavior of the cor- 
relation exponent, and analyse data on the 2 leg t J model 
at electron density n = 3/4, showing that a CDW may 
occur for weak enough J. In Sec. IV we review Friedel 
oscillations in Luttinger liquids and apply this analysis 
to data on the single chain Hubbard and t — J model and 
the 2 leg t — J model. 



II. THE 2-LEG LADDER 

A convenient starting point for bosonization is the 
weak cou pling Hubbard model version of the 2-leg lad- 
der, Eq. (L2) with U/t small. It is then far from obvi- 
ous that this analysis will apply to the infinite c oup ling 
limit, corresponding to the t — J model of Eq. (1.2) so 
that comparisons with numerical results is important. In 
the weak coupling limit we may start by diagonalizing 
the non-interacting problem, giving symmetric and anti- 
symmetric electron operators ipxa where A = e, o labels 
even and odd channels. We pass to the continuum limit 
by introducing left and right moving fields: 



V'LAa(^)+e''"^"^fl.Aa(x) (2.1) 



with kpe and kpo the fermi wave vectors for the two 
bands. In the usual way, we represent the left and right 
moving fermion fields by left and right moving boson 
fields: 



and then, introducing the dual canonical Bose fields, 

(t>\a — 4>R\a + (j)L\a Q\a ~ (t>R\a — (t>L\a (2.3) 

It is then convenient to introduce spin and charge bosons, 
for each channel: 

^Xp = ((^AT +0Ai)/\/2 

(bxa = (<^AT - 0Ai)/\/2 (2.4) 

and then, finally, switch to the two linear combinations 
of the even and odd bosons: 



0±p 



(2.5) 



with similar relations for 4>±a, 0±p, and 6±a- Actually, 
this last transformation is not canonical when the even 



and odd bosons have different isdocities. However, we 
will follow the standard practiccaB of assuming that this 
velocity difference is irrelevant. 

A renormalization group analysis, based on the weak 
coupling Hubbard model, suggests that the cosine inter- 
actions "pin" the bosons 6±^ and introducing ex- 
citation energy gaps for these bosons and leaving c/j-^-p 
as the only massless boson which thus describes the low 
energy excitations. All interactions involving (p+p, O-^-p 
are irrelelvant in this phase. The low energy effective 
Hamiltonian can be written: 



H 



dx 



K+pUlp 



1 



dx 



(2.6) 

where H+p is the momentum density variable canonically 
conjugate to w+p is the velocity of the corresponding 
gapless low energy excitations and the parameter if+p 
controls the correlation exponents.- (Our definition of 
K^p corcfisponds to that of SchulzB and Hayward and 
PoilblanoJ but is the inverse of the_parameter with the 
same name in Balents and Fisher.EI) The Hamiltonian 
may be equally well written in terms of the other boson 
field (j)+p and its conjugate momentum using: 



He — d(j)/dx, = —dO/dx. 



(2.7) 



The two parameters, w+p and i^+p are generally difficult 
to calculate analytically and are extracted from numeri- 
cal data. The long-range behaviors of various correlation 
functions are calculated straightforwardly by expressing 
the corresponding fermionic operators in terms of the 
bosons (/'±CT7 4'±p and their duals. Exponentials of the 
pinned bosons can be replaced by their groundstate ex- 
pectation values but exponentials of the duals of pinned 
bosons lead to exponentially decaying factors in correla- 
tion functions. Exponentials of the gapless boson, 0+p 
and its dual give power law decaying factors. For in- 
stance, to calculate the uniform part of the pair correla- 
tion function we bosonize the pair operator: 



oc e 



i\/7r(0+p-(-0- 



(2.8) 



We may replace the exponentials of 0_p and 9±„ by 
a constant factor leaving simply the operator e*^'^+''. 
The correlation function for this operator decays as 
l/lxj^/'-^^+f''. The same result is obtained for the cor- 
relation function of a pair of electrons in the odd chan- 
nel, Aq. On the other hand, the correlation function 

< Aj(a:)Ao(2/) > is the same except for a factor of 

< e'^^^'I'-p >. This is expected to be < 0, based on 
the sign of the cosine interaction which pins the 4>-p bo- 
son to •\/7r/4, so this pair correlation function has the 
opposite sign, corresponding to "d-wave pairing". The 
2kp part of the ( even, spin-up) density operator is: 



(2.9) 
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The correlation function of this operator has exponential 
decay due to the eT'^^^-p factor. Exponential decay is 
also obtained for all other terms in the Ikp part of the 
density operator. On the other hand, if we consider the 
correlation function for the square of the density opera- 
tor, we get power-law decay for the 4fci? part. This arises 
from terms of the form: 



-2i(fc_Fe+fcFo)a;„/,t 



g-2i(fcFe-|-fcFo)a:g-2iv^(e+p+e+„)^ (2.10) 

We may replace exp(— 2zy^^?+ct) by its expectation value 
leaving only the gapless field which gives power law 
decay: 

< n{.)M^f >oc coB[2(fc^e + fc^o)a: + a] ^ 

In fact we expect that the density operator itself will pick 
up a 4fci? term proportional to this and also exhibit the 
same type of ^kp power law decay. This phase with one 
gapless charge mode and no gapless spin modes, labeled 
"CISO" is expected to be the generic phase in the 2 leg 
t — J or Hubbard ladder. Whether this phase is better 
thought of as almost a d-wave superconductor, or almost 
an incommensurate CDW depends on the magnitude of 
K^p. The pairing correlation function decays less rapidly 
than the density correlation function in the case K^p > 
1/2. Schulz has showrO that K+p ^ 1 as the density, 
n, approaches 1. At n = 1, the ladder is a spin gapped 
Mott insulator. 

Thus it is of interest to calculate K+p. This can be 
done from the finite size spectrum using either periodic 
or open boundary conditions. The latter are often used 
with DMRG so we discuss both cases. While, as we dis- 
cuss below, it is enough to measure the excitation en- 
ergies of only 3 states to extract K+p (and the velocity 
v^p) we give formulas for an infinite number of low lying 
excitations. This allows for a more extensive check of the 
bosonization results with DMRG. Due to the somewhat 
conjectural nature of the extrapolation of weak coupling 
Hubbard results to the t—J model this would be a worth- 
while check. We discuss only the groundstate and low- 
lying excitations, whose excitation energies all scale to 
as 1/L, where L is the length of the system. These can 
be si mply calculated from the free boson Hamiltonian of 
Eq. (|2.6| ) when proper account is taken of the bound- 
ary conditions on the boson field, (?+p. We will generally 
set v+p = 1, restoring it by dimensional analysis when 
needed. 

We first consider the case of periodic boundary con- 
ditions (b.c.'s). Clearly periodic b.c.'s on the fermions 
translate into periodic b.c.'s on the boson fields but we 
must take into account the fact that the boson fields are 
actually phase fields (only their exponentials and deriva- 
tives occur as physical local operators). Since all physical 



operators in the low energy effective theory involve inte- 
ger p owers of ex p\iy/Tr(t)+p\ and exp[i2Y^0+p], as in Eqs. 
.81) and (2.10), we see that: 



9+p{L) = 9iO) + V^p, 



(2.12) 



where m and p are arbitrary integers. To obtain the 
complete low energy spectrum we simply write a mode 
expansion for the boson fields, consistent with these b.c.'s 
and the canonical commutation relations. This gives: 



(f>+pix) =(f)Q + 



2^/TTmx 



E l^ f i2-iTkx/L I -i2Trkx/L , , \ 



fc=i 



ix) 



^/npx 



(2.13) 



k=l 



ink 



(ai^fee''2-fc-/^ - aifeg-^^.''"/^ -f /i.e.) 



Here the qu ant um numbers m and p are integer valued. 
Using Eq. (2/7) we see that the zero wave- vector mode 
of the conjugate momentum, He, has eigenvalues l^jT^m. 
Thus the wave-functions for the zero-momentum modes 
are exp[i2-y/7rm0o]- These wave- functions are invariant 
under ^ + V^j ^ required by the angular nature of 
9. (The same result holds for the dual zero wave vector 
field, 00 its conjugate momentum.) The aj^x and 
o-Bk operators annihilate left and right moving harmonic 
boson modes. We may give a physical interpretation to 
the quantum number m using the bosonization formula 
for the total charge density: 



A,Q 



A. 



^ {-l/V^)Y,d9x,c/dx 

= {-2/^)d9+p/dx. (2.14) 

Thus the total charge, relative to that of the groundstate, 
is: 



Q = (-2/a/^) I dxd9+p/dx = -2A9+p 
Ja 



-2p. 



(2.15) 

Only excitations of even charge occur in the CISO phase. 
The other quantum number, m, measures the "chiral 
charge" which has a less obvious physical interpretation 
as the difference of charges of left and right movers. 

Inserting Eq. (2.14) into the Hamiltonian of Eq. ( |2.6| ) 
and using the identity, Eq. (2.7) we can immediately 
read off the spectrum: 
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E-En = 



-2pii H — ^- 



K+pUi 



AK 



(2.16) 



Here n^fe and njik are the occupation numbers for the 
left and right moving states of momentum ±27rfc/_L. 
is the groundstate energy for a given density, n, and is 
non-universal. This formula gives the excitation energy 
for low- lying excitations with all quantum numbers ^ L. 
The parameters (chemical potential) w+p and K^p all 
depend on density. 

The parameters -f^+p and u+p can be determined by 
measuring the excitation energies of three states. Gen- 
erally one calculates the excitation energies of the states 
with — A(5/2 = p = ±1 (and all other quantum numbers 
set to zero) to determine the ratio v+p/ Kj^p^ using: 



E{p = 1) + E{p = -1) - 2E^ 







■j+p/{K+pL). (2.17) 



The compressibility for a 2-leg ladder is generally defined 
as 



1 d?E 
2L~dn?' 



(2.18) 



where n is the density. From Eqs. ( 2.15 ) and ( 2.16 ), 
this has the value nv+p/ {2K+p). The velocity may be 
measured separately from the excitation energy of the 
lowest state of momentum 2-k/ L: 



E{nn,^l)^Eo^^^. 



(2.19) 



Thus measuring three excitation energies, for fixed, large 
L, allows a determination of the correlation exponents. 
Using a large L is important because corrections to Eq. 
( ^.16 ) are only down by additional powers of 1/L. These 
corrections become especially large at commensurate fill- 
ing near a transition to a CDW. The corresponding expo- 
nent will be given in Sec. III. While these three measure- 
ments are enough to determined the critical exponents, 
measurements of energies of additional states and a study 
of the L-dependence provide additional confirmation of 
the predictions of the RG and bosonization. 

As mentioned above, DMRG works much more effi- 
ciently with open boundary conditions. We now discuss 
the finite size spectrum in that case. If we number the 
sites from 1 to L then we have a free boundary condition 
at J = 1 and j — L. A free boundary condition on the 
fermion field at j = 1 is equivalent to a vanishing bound- 
ary condition at j = 0. [This follows from adding an 
extra "phantom site" at j — with associated hopping 
and exchange terms but then making these terms van- 
ish by imposing ipxa{0) = 0.] In terms of left and right 
movers this free boundary condition becomes: 



1pL\a{0) = -V'-RAa(O). 



(2.20) 



Upon bosonizing we obtain: 



Oxa{0) = constant. 



(2.21) 



(2.22) 



(A determination of this constant involves consideration 
of some subtle commutators. We will not bother to keep 
track of it in what follows.) Clearly this boundary con- 
dition remains the same when expressed in the eventual 
basis p/cr, +/— ■ The only important boundary condition 
for the low energy excitations is: 



0+p{0) = constant. 



(2.23) 



The same b.c. is obtained atx — L + lf^iL except 
that the constant will generally be different. Taking into 
account the periodic nature of 9, Eq. (2.12), and the 



fact that Hg (X d9/dt (X dtp/dx must also vanish at the 
boundaries, we see that the mode expansions become: 



1 



r-^^ \l K+pirk 



■ cos 



fc=l 



Trkx\ ( f\ 



9+p(x) =00 + 



\pii{p — a/27r)x 



L 



+p . . / Trfcx 
/ —I sm 

fe=i ^ 



^ (afe - 4) . (2.24) 



Here a is proportional to the difference of constants ap- 
pearing in the b.c.'s at x = L and x = 0. Substituting 
into the Hamiltonian, we now obtain the finite size spec- 
trum: 



Eq = -2p^l 



nv-i 



L 



{p~a/2TTf 
2K+0 



k=l 



(2.25) 



We see that Eq. ( 2.17 ) for w+p/if+p remains true with 
free b.c.'s as we might expect due to the relation with 
the compressibility. We may now determine the velocity 
independently (and hence determine K^p) by measuring 
the gap to the first excited state with the same charge as 
the groundstate: 



E{ni = 1) - So = 



L 



(2.26) 



Note that this gap has half the value of the gap to the 
lowest energy state of momentum 271 /L in the case of 
periodic b.c.'s, Eq. (2.19). This is just a consequence of 
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the famihar resuh that the spacing of wave-vectors for 
open b.c.'s is tt/L, half the spacing for periodic b.c.'s. 

Previous measurements of the parameter -ftT+p, and its 
analogue in other chain and ladder systems, have used 
periodic boundary conditions. Exact diagonalization, or 
"modified Lanczos" methods work efficiently to deter- 
mined the lowest energy state of given quantum num- 
bers. Thus the energy differences of Eqs. ( [2.17j ) and 
( ^.19 ) can be measured from finding the lowest energy 
states with various charges and with momenta and 
27r/L. An alternative approach is to meaure the depen- 
dence of the groundstate energy (for a fixed charge) on 
an applied flux, i.e. a twist in the boundary conditions 
on the fermion fields: 



(2.27) 



This corresponds to putt ing a twist into the boundary 
condition on (/)-)_p in Eq. ( 2.12| ), as can be seen from the 
fact that all the fermion fields {L or R) contain a phase 
factor of exp[iy/7r(j)^p/2]. 



(f)+p{L) - 0+p(O) + 2^{m + 2$). 



(2.28) 



From Eq. ( 2.16 ) we see that the groundstate energy is 
increased by: 



(2.29) 



Hence measuring the flux dependence of the ground- 
state energy determines v+pK^p, while measuring the 
compressibility determines v+p/K+p. Hay ward and 
Poilblancu measured K+p in this way for the 2-leg lad- 
der using a maximum system size xif 2 x 10. Another 
approach was taken by Siller et al.El They modeled the 
2-leg tJ model as a single chain model of bosonic hole 
pairs. The parameters in the effective Hamiltonian for 
hole pairs were determined from DMRG calculations with 
open boundary conditions on system sizes up to 40 x 2 
for systems containing 2 holes or 4 holes only. The re- 
sulting bosonic model was then studied using exact di- 
agonalization (Lanczos) methods and periodic boundary 
conditions. Since the number of effective bosons (of den- 
sity x/2 where x = 1 — n) is small near 1/2-fllling, it was 
possible to diagonalize large systems (up to 220 sites in 
the case of only 2 bosons). 

The approach advocated here, working exclusively 
with open boundarjf— conditions, has been used previ- 
ously for spin chaingl^ but not, as far we we know, for 
tJ ladders. It allows us, using DMRG, to study much 
larger system sizes, up to 192 x 2, than are accessible 
with Lanczos in the usual formulation.Q It furthermore 
avoids making any assumptions which are necessary in 
treating the hole pairs as bosons, such as the particu- 
lar form of the boson-boson interaction and the absence 
of three-boson terms, which could induce some density 
dependence in the interaction. However, the DMRG cal- 
culations with the accuracy and system length required 
to determine K^p for these ladder systems were surpris- 
ingly difficult. In order to minimize the effect of cor- 
rection terms to the asymptotic formulas, long systems 



were needed. In these systems, one is probing competing 
pairing and charge fluctuations at very low energies and 
large distances, which is difficult in DMRG. Typically 
more than 2000 states were kept in the calculations and 
more than a dozen sweeps were performed. Some of the 
calculations were more accurate and reliable than others, 
howeve r. In Fig. |^ we show a plot of the left-hand side 
of Eq. ( 2.17 ) versus 1/L for three values of (J, n), show- 
ing linear behavior and allowing us to extract v+p/K+p. 
These calculations require only the ground state energy, 
for which an extrapolation in the energy versus the trun- 
cation error is very reliable. Thus this determination of 
Vj^pjKj^p is the most reliable and accurate ingredient in 
the determination of K+p, and, in fact, the results shown 
could be improved with modest additonal effort. 

In Fig. Ijwe plot the left-hand side of Eq. ( 2.26| ) versus 
1/L, again obtaining linear behavior. These calculations 
were difficult. It is necessary to target two states simulta- 
neously, and extrapolation is not very useful. Note that, 
if L is not large enough and the (infinite L) spin gap is 
not small enough, Ag < 2iTv+p/ L, then the first excited 
state with the same quantum numbers might actually 
have = 1, S"^ = and correspond to a spin excita- 
tion, rather than the desired neutral excitation of the -\-p 
field. (It might also be a neutral excitation of one of the 
other gapped boson fields.) In order to check the former 
possibility, we have calculated the spin gap, i.e. the ex- 
citation energy for the lowest state with S*^ = 1 . We find 
that the very nonlinear behavior for J — 0.35, n ~ 0.75 
shown in Fig. |^ is due to a small spin gap in this system, 
of order 0.04. Thus it was necessary to study a 96 x 2 
system, for which we kept 4000 states, and performed a 
dozen iterations, to clearly see the required gap. We be- 



0.25 




0.04 



FIG. 1: These results were obtained by targetting three dif- 
ferent ground states with varying numbers of particles. The 
straight lines are linear fits to the data. The slopes of the 
fits (which asymptotically yield nv/K) are 2.70 for J=.5, 
n=0.875, 6.75 for J=.5, n=0.75 and 7.68 for J=.35, n=0.75. 
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0.00 
0.00 



0.08 



FIG. 2: Excitation energy for lowest state with = and 
same electron density. The spin gap is about 0.14 for J=0.5, 
n=0.875, about 0.10 for J=.5, n=0.75 and about 0.07 for 
J=.35, n=0.75. The breakdown of the asymptotic linear be- 
havior is clearly visible in the third case, at energies of order 
the spin gap. The straight lines are linear fits to the data. 
The slopes of the fits (which asymptotically yield nv) are 
1.631 for n=0.875, J=0.5; 3.03 for n=0.75, J=0.5; and 1.6368 
for n=0.75, J=0.35. 



(n, J) 


V 


-^energy 




(0.875,0.5) 


0.519 


0.604 


0.633 


(0.75,0.5) 


0.964 


0.449 


0.359 


(0.75,0.35) 


0.796 


0.33 


0.284 



TABLE I: Results for the charge velocity v and Luttinger 
liqui d par ameter K. K energy comes from using Eqs. (2.17) 



and (2.26), while Kamp is determined using the decay of the 
Friedel oscilations in the center of the system as a function of 
the system length. 



lieve this data is reliable, but due to the large numerical 
work required we have only studied three different values 
of (J, n). Combining the results of Figs. |l| and |^ allows 
a determination of u+p and hence -fC+p. The resulting 
values are given in Table 0. 



III. CHARGE DENSITY WAVES 

While a CI SO phase is expected over most of the 
parameter range in the 2-leg t — J model, a com- 
pletely gapped, charge density wave (CDW) phase may 
occur at special commensurate filling factors, depend- 
ing on the VjaluCj of J/t (and other possible interaction 
parameters). IIjIIj In this section we review the conditions 
on the critical exponent parameter, if+p for a CDW to 
occur and discuss numerical work on the t — J model at 
electron densities n = 3/4 and n = 1/2, corresponding 
to fillings of 3/8 and 1/4 respectively. 



The low energy effective Hamiltonian only contains 
Fourier modes of the electron fields within a small mo- 
mentum range, ±A of ±kp. Umklapp type interaction 
terms, which do not conserve separately the number 
of left and right movers, are generally accompanied by 
rapidly oscillating phase factors. They therefore do not 
appear in the low energy effective Hamiltonian since the 
electron fields vary slowly and hence these rapidly oscil- 
lating factors cause these interactions to average to zero. 
However, at special filling factors, corresponding to spe- 
cial values of the Fermi wave-vector, the oscillating fac- 
tors become constants and these operators then appear in 
the effective Hamiltonian. Whether or not they produce 
a gap and a CDW depends on whether or not they are 
relevant operators in the RG sense. In a 1-dimensional 
relativistic quantum field theory an operator is relevant if 
it has a scaling dimension, x < 2. This scaling dimension 
determines the exponent, rj with which the correlation 
function of this operator decays, with: 



2x. 



(3.1) 



Using the free boson Hamiltonian of Eq. (p^), the scal- 
ing dimension of any operator is easily obtained. For 
example, the Hubbard or tJ interaction leads to an in- 
teraction term of the form given in Eq. ( 2.101 ). We see 
that this oscillates at wave- vector 2{kFe + kpo)- Only 
when this wave-vector is a mutiple of 27r will this op- 
erator appear in the low energy effective Hamiltonian. 
Furthermore, any operators whose correlation functions 
decay exponentially are irrelevant. 

The number of electrons, N\ in each band is: 



Nx/L = 2kFx/^. 
(The factor of 2 arises from spin.) Thus: 
2{kpe + kpo) — 27rn, 

where: 

n={N,+ No)/2L, 



(3.2) 
(3.3) 
(3.4) 



is the electron density. While the value of kpe and kpo 
may be renormalized by interactions (and may, in fact, 
not really be well-defined in the interacting model) their 
sum is "protected!!, by the one-dim ensi onal version of Lut- 
tinger's theoremllj so that Eq. (3^) is expected to be 
exact. 



Thus we see that the Umklapp operator of Eq. ( 2.1C ) 
can only occur in i?e// for n — I (half-filling). This 
operator has scaling dimension: 



K 



+p- 



(3.5) 



At 1/4 filling, corresponding to n — 1/2, the product of 
this operator times itself with spin up replaced by spin 
down, can occur, giving rise to an operator containing 
8 electron operators. Under bosonization this is propor- 
tional to exjp[—4:iy/n9-f-p\ and has dimension x — AK^p. 



At a filling of 3/8 (or 1/8) corresponding to n=3/4 (or 
1/4) the f ourth power of the basic Umklapp operator in 
Eq. (2.10) containing 16 electron operators can occur. It 
has dimension x ~ IGKj^p. 

The condition on the RG scaling dimension, x, for an 
operator to induce a CDW is somewhat subtle, and de- 
pends on whether the density or J/t is varied. If we 
vary J/t, holding the density fixed at the commensurate 
value, then the system will remain in the gapless CISO 
phase as long as a; > 2, so that the operator is irrele- 
vant. A transition to a gapped CDW phase occurs at the 
value of J/t where x becomes < 2. In the CDW phase 
the value of K^p becomes undefined since all correlation 
functions decay exponentially (or go to constants). We 
may also consider what happens as we vary the density 
in the vicinity of the commensurate value. If the system 
is not in the CDW phase at a commensurate filling, then 
we expect x to vary smoothly and, of course, to satisfy 
X > 2 at the commensurate point. On the other hand, 
if the system is in the CDW phase at a commensurate 
filling, then we expect a; — > 1 as the commensurate filling 
is approached. This, follows from a general and clever 
argument by SchulzEj based on "refermionizing" the sin- 
gle boson Hamiltonian. The relevant operator always 
behaves like a fermion mass term, of scaling dimension 
X = 1 very close to the commensurate filling, whenever 
a CDW occurs at that point. It follows that, x must 
vary rapidly in the vicinity of the critical value of J/t 
where the CDW transition occurs near a commensurate 
filling. On the CDW side of the transition x will be close 
to 2 away from the commensurate filling but then drop 
abruptly to 1 as that filling is approached. We emphasize 
that Schulz's argument, which was formulated in terms 
of the single boson sine-Gordon theory, is very general 
and doesn't depend on the underlying microscopic model. 
Similar behavior can occur for a t — J or Hubbard ladder 
with any number of legs. Such rapid variation of I'C in the 
vicinity of commensurate filling for a system with a CDW 
was observed by Schulz for the single leg Hubbard model 
near half-filling, by obtaining K from the Bethe ansatz 
and similar behavior was also observed near 1/2-filling 
(n = 1) for the 2-leg t — J model by Siller et al. 

Now consider a density n = 1/2 corresponding to 1/4 
filling. If we sit at this density and vary J/t then a 
CDW transition occurs when -fC+p = 1/2, with the gap- 
less phase having if+p > 1/2. On the other hand, if J/t 
is such that the system is in a CDW phase at n = 1/2, 
then K+p ^ 1/4 as n 1/2. 

Similarly, a CDW at 3/8 filling, n = 3/4, is signalled 
by K+p 1/8, as we vary J/t at fixed n or K+p 1/16 
as we vary the density at fixed J/t in the CDW phase. 
See Fig. |. 

The connection between the relevance of the multiple 
Umklapp interaction and the presence of a CDW is eas- 
ily established. When the multiple Umklapp term (for 
n = 1/ 2 or 3/4) i s rel evant it pins the 9+p boson. From 
Eqs. ( 2.10| ) and ( |3.3| ) we see that the 4/cf term in the 
density operator then has a non-zero groundstate expec- 
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FIG. 3: The qualitative behavior of K+p as a function of J/t 
and n. A CDW occurs along the dashed line at n = .75 so 
K+p is not well-defined there. As this line is approached by 
varying n, K+p — > 1/16. On the other hand, at n = .75 and 
larger J/t there is no CDW and K+p is well defined, with a 
value approaching 1/8 at the CDW critical point. 



tation value, so that 

< Tij >oc cos(27™j -|- a) -|- constant, (3.6) 

where a is a constant. Thus the density oscillations 
should have wavelength 2 for n=l/2 and wavelength 4 
for n=3/4. In general, it is sometimes convenient to in- 
troduce 



S = 1 — n, 



(3.7) 



the density of holes measured from half-filling, sometimes 
called the density of "holons" and proportional to the 
doping parameter in the cuprates. We see that n can be 
replaced by 5 in Eq. ( |3.6[ ) so that the wavelength of the 
CDW oscillations is l/S. For a 2-leg ladder, this is the 
average horizontal separation of hole pairs. ( "Horizontal" 
refers to the direction along the legs of the ladder.) 

The CDW corresponds to a broken translational sym- 
metry, so that there are two different groundstates for 
n=l/2 and four different groundstates for n=3/4 (dif- 
fering only by translation by 1 site). For a finite lad- 
der with periodic b.c.'s we expect quantum tunnelling 
between these groundstates to occur so that no density 
oscillations exist in the finite system groundstate. On 
the other hand, with open b.c.'s one or the other of the 
groundstates gets picked out by the boundary conditions 
so that the CDW is directly observable. 

We now return to the point raised in Sec. II, the 
higher order corrections i n 1/ L to the finite size en- 
ergy gaps given in Eqs. ( 2.16| ) and ( |2.25 ). The lead- 
ing correction is determined by the leading irrelevant op- 
erator. At a commensurate filling, when the system is 
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close to having a CDW, the multiple Umklapp operator 
discussed above has a dimension x only slightly greater 
than 2, corresponding to being barely irrelevant. In this 
case the higher order terms in the energy formula are 
of 0{1/ L^~^). This will make it difficult to determined 
K+p reliably from finite size gaps in the vicinity of a 
CDW, complicating the determination of the CDW phase 
boundary. 

Previously reported worki based on the bosonic hole 
pairs approximation together with DMRG and Lanczos, 
on the 2-leg tJ model, with J /t — .35, found a smooth 
behavior of K+p in the vicinity of n — 3/4 with a value 
at a commensurate density n — 3/4 of approximately 
.232. Using the direct DRMG approach, we have found 
K^p « .33. Both results are far above the critical values 
(1/8 and 1/16), indicating that a CDW does not occur 
at n = 3/4 for J/t ~ .35. The charge gap for n — .75 
with J = .35 and J = .25 is plotted versus 1/L in Fig. 
^. There is, indeed, no evidence for a charge gap at the 
larger value of J, consistent with the absence of a CDW. 
On the other hand, there is evidence for a charge gap 
at J = 0.25. Note that these calculations involve only 
ground state energies and thus are very reliable. Al- 
though the results are extrapolated to zero truncation 
error, the extrapolation is quite small in magnitude, of 
order the symbol size in the figure, and the estimated er- 
ror in the extrapolation is very small, as shown by the er- 
ror bars. These results provide some evidence that there 
is a criticial value of J between 0.25 and 0.35 such that 
a charge gap and CDW occur for smaller J. 

As reviewed in the next section, if there is no CDW, 
and the system in is the CISO phase, we expect the den- 
sity oscillations at the center of the chain to decay with 
chain length as L^^+p with X+p > 1/8. Comparing the 
oscillation amplitude for n = .75, J — .25 and the two 
longest lengths studied, L=96 and L=192, we find a ra- 
tio of oscillation amplitudes of .915 which would give an 
exponent of K+p = .128. This is slightly larger than the 
critical value of 1/8. On the other hand, it is so close 
to the critical value that we might wonder if the density 
oscillation amplitude would eventually approach a con- 
stant with still longer chain lengths corresponding to a 
true CDW. 

Another important question is the behavior of the spin 
gap in the putative CDW phase. Since the Umklapp term 
which drives the CDW does not contain the spin bosons, 
it is natural to assume that the gapping of the charge 
boson, corresponding to the CDW, occurs without any 
major effects on the other bosons which could thus re- 
main gapped. However, it is possible that one (or more) 
of the other bosons, such as a spin boson, becomes gap- 
less at the same transition. Indeed, the occurence of the 
CDW may be related to the close proximity of a phase 
with no spin gap. 

In Fig. |[ we show the spin gap on a 16 x 2 system for 
a wide range of J. The data suggests that in the thermo- 
dynamic limit, the spin gap vanishes near J = 0.25. To 
verify this result, one must perform a finite size study. 



0.6 




FIG. 4: Scaling of the charge gap with 1/L, showing the 
occurance of a CDW at n = .75 for J — .25 but not J — .35. 



However, the spin density pattern on the 16 x 2 system 
for J — 0.2 has the largest values of (5f) concentrated 
near the ends of the system, indicating that an edge exci- 
tation has the lowest energy. In order to ensure that we 
have a bulk spin excitation in the calculation, we next 
studied systems where we increased the rung exchange 
interaction to J + 0.3 on the first and last rungs of a lad- 
der. This drives up the energy of the edge excitation, and 
one finds that the resulting spin pattern is concentrated 
in the central region of the ladder. In Fig. ^, we show the 
spin gap as a function of 1/L for several systems. The 
results for J = 0.2, n ~ 0.75 are consistent with a gapless 




FIG. 5: Spin gap on a 16 x 2 system as a function of J/t. 
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FIG. 7: A charge arrangement with period 2 at n=3/4, cor- 
responding to a gapless CDW. The black circles represent 
electrons and the white circles holes. The ovals indicate that 
the location of the hole is symmetrized between the two sites 
on a rung. The solid lines indicate the formation of dimer 
singlets. The arrows indicate the spins of the unpaired elec- 
trons, but we donnot expect them to be Neel ordered in this 
one-dimensional system. 
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FIG. 6: Spin gap as a function of l/L for several systems. 
The values of the exchange on the first and last rungs have 
been altered to avoid edge excitation. 



FIG. 8: A charge arrangement with period 4 at n=3/4, corre- 
sponding to a gapped CDW. (The actual density oscillations 
seem to correspond to pairs of holes shared evenly by every 
second pair of rungs.) 



state, although one can never rule out a very small gap. 
Similarly, for J = 0.25, n = 0.75 no evidence for a spin 
gap on the largest system, a 192 x 2 ladder. For J = 0.2, 
n = 0.875 we see a clear spin gap of order 0.013. 

The proximity to a phase with vanishing spin gap may 
be related to the occurence of the CDW. As the spin gap 
gets smaller, and the corresponding length scale larger, 
we may expect the size of the 2-hole pairs to get larger. 
This is naturally associated with a growing scattering 
length for the effective boson system, which enhances the 
formation of a CDW. 

A CDW with vanishing spin gap dX n — 3/4, can be 
easily understood heuristically. One can imagine one hole 
localized on every second rung, with the electrons on dou- 
bly occupied rungs forming spin singlets, as illustrated in 
Fig. 1^, giving an effective 5 = 1/2 Heisenberg chain with 
lattice spacing 2 which has vanishing spin gap. On the 
other hand, in the spin-gapped phase we can think of 
pairs of holes localized on every fourth rung, as in fig. ||. 
Thus we expect the groundstate density oscillations to 
have wavelength 2 in the CDW with vanishing spin gap 
but wavelength 4 in the gapped CDW. 

The connection between wavelength doubling and t]*© 
spin gap is related to the Lieb-Schultz-Mattis theorem.Ej 
If the wavelength is not doubled to 4, then gapless exci- 
tations are expected. More precisely, we can prove the 
existence of gapless excitations at n=3 /4 if the ground- 
state is invariant under translations by 2 and under site 
parity: j — > —j and is not ferromagnetic. (Note that 
link parity, j — > —j + 1 is spontaneously broken by the 
wavelength 2 groundstate illustrated in Fig. ^ but site 
parity is not.) This is proven by considering a long lad- 
der of even length L with periodic b.c.'s. Let > be a 
groundstate. We then consider the variational state ob- 



tained by acting on the ji/'o > with the unitary operator: 
C/ = exp[z(2VL)^jn,T], (3.8) 



where riji is the number operator for spin up electrons 
on rung j (summing over both sites on the rung). It is 
straightforward to calculate < ■0o|C^^(^f — £'o)IV'o > and 
show that it is 0(1/L) provided that the following term 
vanishes: 



3 



[< Mc]x^Cj+ix^\i}jo > - < V'o|c]+iat'=jmI'/'o >]■ 



(3.9) 

We can show that this vanishes, for example, by using 
the assume site parity symmetry to show that: 

< i^o\c]x^Cj+ix-\\i}o >=< ■(/'o|c];^|Cj_iatI^/'o > • (3.10) 

To complete the proof of a low energy state we must 
prove that U\il>o > does not become \ipo > in the limit 
L ^ oo >. We prove this by considering the behavior of 
U under translations by 2 sites, the assumed symmetry 
of the groundstate. We fine that under this translation: 



(3.11) 



where n| is the density of spin up electrons, n| — iV| /2L, 
where iV| is the total number of spin up electrons and 2L 
is the number sites on a 2-leg ladder of length L. If we 
further assume that rt| = n^, i.e. that the groundstate is 
not ferromagne tic, t hen we may replace the exponential 
factor in Eq. ( |3.1lD by e*''''". For a density, n = 3/4, 
this factor is -1. This proves orthogonality oi U\^o > 
with 1-00 > and hence the existence of a low energy ex- 
citation. On the other hand, if the groundstate is only 
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invariant under translations by 4 sites then, repeating the 
argument with a translation by 4 sites changes the expo- 
nential factor to e*^'^" which is 1 for n=3/4. Thus we 
cannot prove that U\^po > is orthogonal to the ground- 
state in this case; it may approach it in the limit L oo, 
so the proof collapses. Unfortunately, this does not prove 
that there is a gap when the groundstate has wavelength 
4, only that the gap neccessarily vanishes when it has 
wavelength 2. However, it often appears to be the case 
that when the LSM theorem fails, a gap appears, for 
instance, for S=l/2 spin chains, a gap is expected when- 
ever the groundstate has wavelenth 2, as in a dimerized 
state. Thus, the nature of the density oscillations pro- 
vides additional evidence for which phase the system is 
in. 

To see what Fourier components are present in the den- 
sity oscillations, we Fourier transform the density as a 
function of rung position and plot the power spectrum. 
In order to avoid spurious edge effects, a windowing func- 
tion is used. If the original density is d{x), we Fourier 
transform D{x) = W{x)d{x), where W{x) is a smooth 
windowing function. W{x) is chosen to vanish at a; = 1 
and X — L, and to be unity for the middle third of the x 
values, with a smooth continuation in between; the par- 
ticular choice we use is given in Ref . (|l^) . Before Fourier 
transforming, we calculate the sum / of D{x) from 1 to 
L, and the sum J of W{x), and then subtract from D{x) 
the function J VF(a;). This removes a large peak centered 
at A: = which is uninteresting. The lattice spacing is 
set to 1, so that the exponential in the FT is exp(— ?fcj), 
where j runs over integer lattice sites. With this approach 
there is no restriction on the allowed values of k; the finite 
value of L instead leads to finite widths for the peaks. 

The density oscillations for J = .25 and also J = .2, 
with n = .75 have a very strong component at a wave- 
length of 4, corresponding to a spin-gap, as shown in 
figure ^ 



IV. A POSSIBLE PHASE DIAGRAM 

In this section we briefly discuss a possible phase dia- 
gram and the qualitative behavior of -fsT+p as a function 
of J and n for the 2-leg t — J model. We are interested in 
the region rt > 0.5 where both the even and odd parity 
bands have carriers in the weak coupling Hubbard limit. 
At still lower doping the system is expected to enter a 
CI SI phase. In the region n > .5 our phase diagrarn. 
is similar to that proposed by Hayward and Poilblanc,Ll 
and is based on their exact diagonalization results, earlier 
DMRG results, our new DMRG results and some general 
results which follow from bosonization and RG. This is 
sketched in Fig. |l^. Here the two ordered CDW phases 
at n=0.75 and n=0.5 are indicated by solid cuts for small 
values of J/t and a phase separation region occurs for 
large J/t values. We conjecture that the remaining re- 
gion is in a CISO (d-wave-4ki? CDW) phase. As discussed 
in Sec. II, this phase has the 0_p and bosons pinned 
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FIG. 9: Smoothed Fourier transform of the density oscilla- 
tions showing only a large peak at wave- vector 7r/2 corre- 
sponding to a wavelength of 4. 



and is characterized by power law d-wave pairing and 
4ki? CDW correlations. When K^p > 0.5, the pairing 
correlations are dominant. Supporting this conjecture, 
DMRG calculations for J/t and n values near the phase 
separation boundary show clear power law d-wave like 
pairing correlations. 

The phase diagram shown in Fig. |l^for the t—J ladder 
differs in several ways from that expected for the 2-leg re- 
pulsive U Hubbard model. Firstly, it is generally believed 
that the Hubbard ladder doesn't exhibit phase separa- 
tion and the ordered CDW phases are absent. Secondly, 
weak coupling renormalization-group studies, bosoniza- 
tion [4] and DMRG calculations suggest that K+p < 1 
for the repulsive U Hubbard model. For an extended 
Hubbard model with additional interactions |One can have 
K^p > 1; however Orignac and Giamarchc argue that 
when K+p > 1, the d-wave-4ki? CDW sector gets re- 
placed by a CISO s-wave pairing or an orbital antifer- 
romagnetic phase. This raises questions regarding the 
phase diagram show in Fig. |l^. We believe that while 
alternate CISO phases are certainly logical possibilities, 
they do not occur in the t—J ladder. The argument that 
they occur whenever K+p > 1 appears to depend upon 
a weak coupling analysis which may not be applicable to 
the t—J ladder. In this analysis, the phase boundaries 
between the various CISO phases is determined by the 
equality of 2 small marginal coupling constants in the 
Hamiltonian. These coupling constants also determine 
-R"+p, in the weak coupling limit, in such a way that when 
they are equal, K+p = 1. [See Eq. (35) of Ref. (§.] For 
stronger coupling, and, in particular, for the t—J model, 
the possible phase boundary to a phase with orbital anti- 
ferromagnetism need not be related in any general way to 
the value of K+p and the DMRG results suggest that or- 
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FIG. 10: Schematic phase diagram for the t — J ladder as 
a function of J/t and the electron density n in the region n 
greater than about .5. 



bital antiferijamagnetism does not occur in the ordinary 
t - J model.Ell 

The phase separation boundary was obtained in a pre- 
vious DRMG study.E3 As this boundary is approached 
the compressibility diverges and -ftT+p — > oo. 

We have taken into account the behavior of i^+p near 
commensurate fillings in drawing Fig. 
at n ^ 1, the t — J model reduces to the 2-leg Heisen- 



berg spin ladder and has an infinite gap for all charge 
excitations and presumably a gap of 0{ JXJot spin exci- 
tations. The general arguments of SchulzaO then imply 
that K+p — > 1 along the line n = 1. We have also as- 
sumed that a CDW occurs at n—.75 and n=.5 for small 
enough J. As discussed in the previous section, we have 
found some evidence for such a phase at J — .25 and 
n=.75. The corresponding behavior of K^p, reviewed in 
the previous section is incorporated in Fig. |l^. Note 
that, at small J, -fC+p apparently varies rapidly between 
.0625 for n=.75 and 1 for n=l. As we also discussed in 
the previous section, it is possible that the spin gap van- 
ishes at commensurate fillings for some range of J. We 
do not attempt to indicate this possibility in Fig. 10. 

For values of J/t around .35, which has been ar- 
gued to be a reasonable value for modeling the cuprates, 
if+p < .5 near n=.75 but increases towards 1 as n=l 
is approached. Thus there is only a narrow window of 
doping near 1/2-filling where d-wave pairing correlations 
dominate. 



V. FRIEDEL OSCILLATIONS 

At incommensurate filling, or in general when no CDW 
occurs, there will still be density oscillations produced by 
the boundaries of an open ladder which only decay slowly 
(with a power law) into the ladder. Detailed predictions 



can be made about these using bosonization generalizing 
the approach of Ref . (||,||) for the spin chain and Hubbard 
chain to the case of the ladder. These oscillations provide 
a further check on the theory. Thus consider again the 
Akp term in the density operator, of Eq. ( 2.1C| ): 



Hi 



Ae 



-2i-Knj —2i^/TT9 



h.c. 



where A is 



< e 



a constant into which we have 
>. As we discussed in Sec. II, 



(5.1) 

adsorbed 
this Akp 

term in the continuum representation of rij, determines 
the leading behavior of the density-density correlation 
function at long distances since any 2k p terms decay ex- 
ponentially. For an infinite system, the correlation func- 
tion of the exponential operator appearing in Eq. ( ^.l| ), 
behaves as: 



< e 



'{r)e 



-2i^e 



+p(0) c\r\-^^+'' 



(5.2) 



2K+ 



where c is a constant, with dimensions of (length) 
Decomposing 6 = (pB, — (t>Li we may write this correlation 
function as a square of two identical factors, the correla- 
tion functions of exp[2iy/Tr(pfjj]^]. It thus follows that the 
oscillating term in the density correlation function is: 



< n{r)n{0) cos(27rm- + a)2c|ylp|r| 



-2K, 



(5.3) 



Now consider the semi-infinite system (r > 0) with one 
free boundary condition. As remarked in Sec. II, this 
corresponds to a b.c. on the boson field given in Eq. 
In particular, ( p.23D , or equivalently: 



(j)Lit, 0) — (f>R{t, 0) + constant. 



(5.4) 



Now using the fact that 0^ is a function ofvt—r only and 
a function oivt + r only, we see that this equation im- 
plies that we may regard as the analytic continuation 
of to the negative axis: 



4>R{r) — (j)L{—r) -f constant. 



(5.5) 



Thus the expectation value of n{r) reduces essentially to 
the square root of the correlation function < n{r)n{Q) >: 

< Uj >^ 2^\A\ cos(27rnj -I- /3)(2j)~-^+''. (5.6) 

The wave- vector of the Friedel oscillations is the same as 
the wave-vector governing the long-distance oscillations 
in the correlation function. No Friedel oscillations occur 
at wave- vector 2kFi as they would for the non- interacting 
system or more generally in phases where, for example, 
all 4 bosons are gapless. The exponent governing the 
decay of the Friedel oscillations is A'+p, 1/2 the expo- 
nent governing the decay of density correlations. Fur- 
thermore, the amplitude is simply the square root of the 
amplitude of the density correlation function (multi plied 
by 2^/^~^). Finally, we may readily generalize Eq. ( |5.6| ) 
to the case of a finite chain of length L, by a standard 
conformal transformation: 



2^\A\cos(2imj + [}) 
[(2L/7r) sin(7rj7L)]^+p 



(5.7) 
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FIG. 11: The amplitude, A, of the Friedel oscillations in the 
center of the system, as a function of the length, L. The 
resulting slopes are — if+p and are given in Table The 
triangles show results for L = 64 as a function of the number 
of block states kept, m, for J = .35, n = .75. One can see 
the very slow convergence of the amplitude with the number 
of states kept. 



The decay of Friedel oscillations allows a way of deter- 
mining the critical exponent K+p which is alternative 
to measuring directly the long distance behavior of the 
density correlations. Indeed this latter measurement be- 
comes quite difficult with open b.c.'s due to the necessity 
of eliminating boundary effects. 

We emphasize that the wave- vector of these Friedel os- 
cillations 2(kpe + kpo) = 27rn, or equivalently 2'kx where 
a; = 1 — n is the hole density, is a characteristic feature of 
the particular CISO phase that we are assuming. More 
correctly, this is the minimum Friedel oscillation wave- 
vector since higher harmonics are also expected to occur. 
This minimum oscillation wave- vector would be different 
in a different phase. This wave-vector corresponds to 2 
holes per wave-length and is the same wave-vector that 
would occur for a 1-component spinless hard core bose 
gas, which is an approximate description of the CISO 
phase in which holej^pairs are assumed to form tightly 
bound rung singlets.Q 

An immediate consequence of Eq. (5.7) is that the 



amplitude of the density oscillations near the center of 
a finite chain scales as L~^+''. We show a log- log plot 
of this amplitude versus length in Fig. |ll|. Fitting to a 
straight line allows a determination of K+p. The corre- 
sponding values of K^p determined in this way are shown 
in Table ||. We see that they are roughly comparable to 
the values obtained from the finite size spectrum. 

In Figs. |l2| and |l^ we show the Friedel oscillations in 
the density at site /, for two different values of (n, J), 
fitted to Eq. ( pT7| ) with i^+p taken as a free parameter. 
The agreement is fair although the presence of corrections 



due to irrelevant operator effects is evident. In both cases 
the value of if+p so determined is in rough agreement 
with the values in the third column of Table I, determined 
from the finite size spectrum. 

We emphasize that these formulas are true very gen- 
erally for Luttinger liquids with a single gapless charge 
boson. In particular, they apply to the spinless single 
chain model. In the non-interacting case we may readily 
find the exact formula for the Friedel oscillations. For N 
electrons on L sites: 



N 



N- 



L 
1/2 



■m— 1 

sin[27rj(A^ 



2 TTTOJ 



L + 1 

+ l/2)/{L- 



1)] 



L + 1 2(L + l)sin[(7rj7(L + l)] 

For large A'^ and L this can be approximated: 

sin 2Tmj 
' 2L sm{nj/L) 



(5.8) 



(5.9) 



This has the expected form of Eq. (5.7) with K = 1 and 
Y^|yl| = l/27r. On the other hand, the 2kp part of the 
density correlation function at long distances is: 



cos 27rnj 
27r2|j|2 



(5.10) 



which has the form of Eq. ( ) with the same value of 
V^\A\ = 1/2^. 

We see that, not surprisingly, when the charge density 
correlations drop off slowly, so do the Friedel oscillations. 
In particular, this makes it difficult to determine numer- 
ically whether or not a CDW occurs at n = 3/4, for ex- 
ample, by measuring density oscillations. If J/t is such 
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FIG. 12: Density at site / from DMRG (* symbols) compared 
to Eq. (5.7) (circles and lines)for (n, J) — (.5, .875) using 
K+p = .63. 
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FIG. 13: Density at site I from DMRG (* symbols) for 
{n,J) = (.35, .75) compared to Eq. (5.7) (circles and lines) 
using K+p = 0.33. 



that the system almost has a CDW then i^+p will be 
only slightly greater than 1/8. The extremely slow decay 
of the Friedel oscillations will be difficult to distinguish 
from a true CDW, where the oscillation amplitude goes 
to a non-zero constant far from the boundaries. 

We now consider the single chain Hubbard and tJ mod- 
els. The Hubbard model is known to be in a ClSl phase 
for all densities (except n = 1) and all U > 0. Thus the 
low energy effective Hamiltonian contains both spin and 
charge bosons. (We use the same notation as for the 2- 
leg ladder except that the ± labels are no longer needed 
since there is only one type of charge boson and one type 
of spin boson.) The charge boson Hamiltonian is written 
exactly as in Eq. ( |2.6D , in terms of the parameter Kp. 
The spin boson Hamiltonian also has exactly this form 
but with = 1, as follows from SU{2) spin rotation 
invariance. The pairing, 2kF and Akp density operators 
now take the form: 

Note that, in this case, we have formed a Akp operator 
in which the spin boson doesn't appear. It then follows 
that the pair correlations, 2k p density correlations and 
Akp density correlations decay with power law exponents: 



'^pair ^ 

V2kp = 1 



= 4i^„ 



(5.12) 



Open b.c.'s imply: 
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FIG. 14: Smoothed Fourier transform of the density, rij in the 
single chain Hubbard model at 7/8 filling for various values 
of the repulsion strength, U. 
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constant. 



(5.13) 



It thus follows that the density will exhibit both 2kp and 
4k p Friedel oscillations with different exponents: 

^ A cos(7rnj + a) ^ B cos(27mj + /?) 



J 



(l+-ffp)/2 



J 



2K„ 



(5.14) 

Where we have expressed the oscillation wave- vectors in 
terms of the electron density using the free electron result 
for electrons with spin on a single chain: 



2kf 



nn. 



(5.15) 



constant 



The amplitudes, A and B are proportional to the square 
roots of the corresponding terms in the density correla- 
tion function and depend on density and U (as does Kp). 

It is interesting to consider the limit U >> t where the 
no double occupancy constraint is present. If we consider 
some initial configuration of holons and spins, then the 
order of the spins along the chain can never change under 
nearest neighbor hopping processes consistent with no 
double occupancy. It follows that the charge dynamics 
must be identical to those of a system of non-interacting 
spinless fermions with the same density. That is to say, 
the charge dynamics is "spin blind" . It follows that the 
density oscillations should not have the nn term in this 
limit but only the 2Tm term as for spinless fermions [Eq. 
(^.9|) ]; i.e. that ^ ^ at C/ — > oo. This argument also 
determines Kp = 1/2 and B = l/27r in this limit. We 
remark that this argument is special to the single chain 
model with only nearest neighbor hopping. On a ladder, 
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or even on a chain with next nearest neigbor hopping, 
spin rearrangement is possible without double occupancy 
by moving the electrons around or past each other. 

In Fig. 14 we plot the Fourier transform of the density 
for the Hubbard model, at 7/8 filling, for various values 
of U. The same smoothing procedure was used that we 
mentioned in Sec. III. In this case 2kp = nn = Ttt/S 
and Akp = 7tt/A or equivalently 7r/4. Both 2kF and 4kp 
components are clearly visible and it is evident that the 
2kp component vanishes at large U. 

What does this imply about the (single chain) t — J 
model? The J ^ limit of the t — J model is identical 
to the [/ — > oo limit of the Hubbard model. Therefore, 
in that limit, the density oscillations have vanishing 2kp 
component. It is not completely obvious what happens 
at non-zero J. Insofar as this is the same as large but 
finite U, we would expect to have a small 2kp compo- 
nent to the density oscillations. Note however, that the 
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no double occupancy constraint is enforced exactly in 
the t — J model, by construction, whereas it is not in 
the finite U Hubbard model. Thus we must understand 
whether it is the no double occupancy which is respon- 
sible for the vanishing 2kp oscillations or whether the 
vanishing of all spin rearrangement processes is necces- 
sary. The argument in the previous paragraph seems to 
rely on the vanishing of spin rearrangement processes so 
we suspect that the Friedel oscillations should exhibit a 
2kp component for any finite J. Of course, for small 
J/t this is expected to be small. In Fig. 15 we show 
the Fourier transformed density for the single leg t — J 
model at density n = 1/4 for two values of J/t. Note 
that now 2kp = ir/A and Akp = tt/2. For smaller J/t the 
4kp oscillations clearly dominate although a small 2kp 
part is observed. At larger J/t, the oscillations at 2kp 
are larger than those at Akp. 



Thus we expect tJ models to exhibit generic Friedel 
oscillations in all cases except for the special case of a 
single leg in the limit J/t ^ Q where the absence of 
spin rearrangement processes eliminates the lowest wave- 
vector component of the oscillations. 



The Friedel oscillation wave-vector can be a useful di- 
agnostic of which phase a particular multi-leg tJ model 
is in. For instance, the minimum oscillation wave-vector, 
2kpe + 2kpo = 27rn, is a characteristic of the particular 
CISO phase of the 2 leg ladder reviewed in Sec. II (in 
which the 9±a- and (jj^p fields are pinned). Other phases 
exhibit other oscillation wave- vectors. We expect such 
an analysis of DMRG results to be useful in determining 
the phase diagram of multi-leg ladders. 
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FIG. 15: Smoothed Fourier transform of the density, rij in the 
single chain tJ model at 1/4 filling for a) J/t = .35 and b) 
J=1.0, showing Friedel oscillations at 2kF ~ 7r/4 and 4fcF ~ 
n/2. 
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